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ABSTRACT: 

The distribution of the Stirling numbers S(n,k) of the second kind 
with respect to k has been shown by Harper [Ann. Math. Statist., 38 
(1967), 410-414] to be asymptotically normal near the mode. A new single- 
term asymptotic representation of S(n,k), more effective for large k, is 
given here. It is based on Hermite's formula for a divided difference 
and the use of sectional areas normal to the body diagonal of a unit 
hypercube in k-space. A proof is given that the distribution of these 
areas is asymptotically normal. A numerical comparison is made with the 
Harper representation for n=200. 
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1. Introduction. 



Previous asymptotic representations of Stirling numbers S(n,k) of 
the second kind have been of two types. One type has been a complete 
infinite series expansion as given by Hsu [1], and by Bleick and Wang 
[2] and [3]. A second type has been the single-term representation of 
S(n,k) given by Harper [4] as the normal distribution approximation 



( 1 ) 



S(n,k) ^ — — exp [- (k-y) ^/2a^] 

a/Jrr 



where the mean y and the variance a are expressed in terms of the Bell 

numbers B by 
n 



( 2 ) 



y 



n+1 



/B 

n 



1 



and 

(3) 



2 



a 



n+2 



/B 

n 



«WV 



1 . 



The purpose of this note is to give a new single-term asymptotic re- 
presentation based on Hermite T s formula for a divided difference, and to 
compare it with that of Harper. 

2. Use of Hermite f s formula. 

A Stirling number S(n,k) of the second kind is defined as the kth 
difference of z 11 at z=0 divided by kl • By [5,p.l0] we find that this 
divided difference can be represented by a formula of Hermite as the re- 
peated definite integral 

(4) S(n,k) = / dt 1 / 1 dt .../ k 1 (d k u n /du k )dt, 

0 0 0 1 IK 

where u^=t^+t 2 +* • +t j c * We imagine that t^, t^ constitute a set of 
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rectangular Cartesian coordiantes and impose an orthogonal transforma- 
tion of coordinates to u^, u^. The volume of the space over 

which the integration in (4) is performed is a portion of a unit hyper- 
cube in k-space. If we allow the coordinate u^ to vary along the body 
diagonal of the hypercube from 0 at one vertex to k at the opposite 
vertex, the sectional areas normal to the diagonal cut by the hyper- 
plane u 1 =t 1 +t 0 +. .+t from the domain of integration define a positive 
function g(u^,k) even with respect to the argument u^-k/2. We take the 



integral of g(u^,k) to be 

(5) 



/ g(u ,k)du = l/k! 
0 



to agree with the volume of the space over which the integration in (1) 
is performed. We drop the u^ subscript henceforth. Noting that g(u,k)=0 
for k<u<0, we find that 

(6) g(u,l) = 1 for 0 £ u <_ 1 , 

(7) 2!g(u,2) = (1 - | u-1 | ) for 0 £ u £ 2 , 

and 0 

(8) 3! 



Hg(u,3) = | 



(3/ 2- | u-3/2 j ) / 2 for 1/2 £ | u-3/2 1 £ 3/2 
2 



3/4 - (u-3/2) for 1 u 2 

Consideration of the Laplace transforms of (6), (7) and (8) suggests that 

we conjecture the Laplace transform of k!g(u,k) to be 

(9) (l-e -s ) k /s k = e -ks/2 ( sinhs/2 ) k 

s / 2 

for all k. We demonstrate the truth of this conjecture later. On perform- 
ing the integration in (4) over the variables u , u , u we find 

Z. K. 



( 10 ) 



S(n,k) = k! (^) j u n k g(u,k)du . 
k 0 
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Using operation 82 of [6,p.l0] on the Laplace transform of 

u 

(11) k! / u m g(u,k)du 

0 

we find the mth moment of the k!g(u,k) distribution about u=0 to be 



( 12 ) 



lim (-l) m (d/ds) m (l-e S ) k /s k . 
s-K) 



It is now easy to demonstrate the truth of the conjecture (9) by show- 
ing, with the aid of the multinomial theorem, that (12) is the same as 
the repeated integral 



(13) 

o 0 

over the volume of the hypercube 



/ dt i/ dt 2 .../ ( t i + t 2+ .. + t k )'"d S 



0 



, 2 2 , , ks 2 /24 , (ks 2 / 24) 2 
exp (a s / 2) = 1 + — ^7 — + “ — JT + • • 



Use of (12) and (5) shows the variance of the k!g(u,k) distribution 
to be 

(14) o 2 = k/ 12 . 

Using (14) the series 

(15) 

is the bilateral, but not s multiplied, Laplace transform of the 
normal distribution 

(16) (l/a/2iF)exp (-t 2 /2a 2 ) 

according to [7,p.2J]. The corresponding series for (9) multiplied by 
ks / 2 

e , or the bilateral Laplace transform of k!g(u,k) shifted left by 



k/2, is 
(17) 



(2/s) k sinh k s/2 = [ 



. u + + (s 2 /4) 2 k 

1 3! 5! J 



2 n 

The dominant k power term in the coefficient of (s /4) in (15) is 
k n /6 n n!, and may be shown to be the same in the expansion of (17) 
by the use of the recurrence formula 6.361 of [8,p.ll9]. This proves 
that the k!g(u,k) distribution is asymptotically normal as k-*». It is 
remarkable that the normal distribution should arise in the purely 
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geometrical context of sectional areas normal to the body diagonal of 
a hypercube of high dimension. 



On replacing k!g(u,k) in (10) by its Gaussian normal approximation 

2 

of mean y=k/2 and variance a =k/12 we find 



(18) 



oo 

S(n,k) 0/ — — (™) | u n ^exp[-(u-k/2)^/2a^]du 



a/2iT 



0 



i /3k , 2 

„ 1 ,n r . .n-k -t /2 Jj 

^ — (.) / (k/2-at) e dt . 

/2tT -°° 



3. Numerical example. 

Table 1 compares the exact values of S(200,k) with the asymptotic 

approximations computed from the single-term representations (1) and 

2 76 

(18). Harper's representation (1), which uses B^qq=. 62475 10 ’ 



y=49 .975 


and a=3.0551, gives 


an excellent fit near the 


mode (k=50) , 


but (18) 


gives a much better 


fit for large values of k. 






Table 1. 


Values of S(200,k) 






Asymptotic 




Asymptotic 


k 


from (1) 


Exact 


from (18) 


2 


222 

.23135 10 


.80347 10 6 


.69244 10 126 


40 


273 

.39504 10 Z/J 


273 

.24458 10 


273 

.42658 10 


50 


275 

.81579 10 


275 

.81493 10 


277 

.15285 10* 


60 


273 

.37452 10 


273 

.53533 10 


274 

.29658 10 


100 


217 

.49065 10 


235 

.22839 10 


235 

.27994 10 


150 


43 

.13938 10 


143 

.30251 10 


.30441 10 143 


199 


.16955 10 -241 


.19900 10 5 


.19900 10 5 
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